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• Abstract 

' Pseudo-Hermitian operators can be used in modeling electromagnetic wave propagation 

in stationary lossless media. We extend this method to a class of non-dispersive anisotropic 
Q^' media that may display loss or gain. We explore three concrete models to demonstrate the 

utility of our general results and reveal the physical meaning of pseudo-Hermiticity and quasi- 
Hermiticity of the relevant wave operator. In particular, we consider a uniaxial model where 



^ I this operator is not diagonalizable. This implies left-handedness of the medium in the sense 

that only clockwise circularly polarized plane-wave solutions are bounded functions of time. 
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1 Introduction 

^ ■ Consider a stationary (time-independent) dispersionless and source-free medium with permittivity 
■ and permeability tensors £ = £{x) and fl = fl{x). Then the Maxwell's equations are given by [1]: 

V ■ (eE) = V ■ B = 0, (1) 
E = £ D/i B, (2) 
B = (3) 

where an overdot stands for a time-derivative, and T) is the curl operator, e.g., DE := V x E. As 
shown in Ref. [2], Eqs. ([2]) and (E]) are equivalent to 

B{f, t) = Mr) - [ ^E{r, r) dr, (4) 
Jo 

E + n^E = 0, (5) 
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where -Bo(^) := -8(^^,0) and 

Vl^:=E ^11 D. (6) 

We can view 2) as a linear operator acting in the Hilbert space H of square-integrable complex 
vector fields F : M'^ — t- C. This is defined by the L^-inner product, 

-< F\G y:= I F{xy ■ G{x) 

where a dot placed between F* and G stands for the usual dot product 

It is easy to see that 2) is a Hermitian operator acting in T-L. Furthermore, for the cases that £ 
and H are given by invertible Hermitian matrices, they also define invertible Hermitian operators 
acting in "H. This in turn implies that Vt^ satisfies the pseudo-Hermiticity relation [3]: 

In particular, for a lossless (and gainless) medium £ is a positive-definite operator, and Vl? is 
quasi-Hermitian i.e., it is related to a Hermitian operator via a similarity transformation [S]. 
This observation suggests an interesting spectral method of solving the wave equation ([S]). Ref. [2] 
gives a rather explicit application of this method for an effectively one-dimensional model involving 
localized inhomogeneity. Ref. [6] studies the possibility of extending this method to a class of 
lossless dispersive media. 

The purpose of the present article is to use the spectral properties of the operator in modeling 
electromagnetic wave propagation in stationary non- dispersive anisotropic media that display loss 
or gain. This requires considering non- Hermitian permittivity and permeability tensors [3 [8] . In 
this case, is generally no longer quasi-Hermitian, and a direct application of the spectral method 
outlined in Ref. [2] is intractable. In the following we will examine the structure of Vt^ for a lossy 
anisotropic medium. Our investigation is motivated by the well-known fact that in many physical 
applications one deals with such materials. A typical example is a magnetic photonic crystals 



2 Solution of the Wave Equation 

The wave equation ([5]) admits the following formal solution [2] 

E{r, t) = cos{Qt)Eo{f) + Q'^ sm{Qt)Eo{f), (8) 

— * — * — # — # 

where Eo[f) := E{r,0), Eq{x) := E{r,0), and 

In this section we will study the structure of the solution ([8]) for the cases that the permittivity and 
permeability tensors, £ and /i, are given by constant non-Hermitian 3x3 matrices. We achieve 
this by performing a Fourier transform. 

^For allFe'H , -<F|_F;^<oo. This in particular implies that |F(a;)| -> as |x| — !• od. 
^See [To] for a comprehensive discussion of lossy crystals. 
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In the Fourier basis the curl operator D is represented by the antisymmetric Hermitian matrix: 




2)= -ih iki \ , ki,k2,heR, (10) 



while £ and /i that appear in the expression ([6]) for f2 are 3x3 complex matriceslil It is 
not difficult to compute in the Fourier basis. This gives rise to a 3 x 3 matrix involving kt and 
the entries of £ and /i, i.e., a total of 39 real parameters. Rather than giving this complicated 
expression, we will examine the Jordan canonical form of Q"^. In fact, it is more convenient to work 
with the dimensionless operator. 



■.= ujQ^n^ = ujq^£ ^D, (11) 



where 



Uo:=c\k\ = cJkl + kl + k1, k := \ k2 . (12) 




In terms of (ED and (191) read 



E{f,t) = cos{QuJot)Eo{f) + Uq^Q-^ sm{nuJot)Eo{r), (13) 
cos(fi uot) := J2 Ht'^'T. (14) 



n=0 



n=0 



A straightforward consequence of (|TOl) is = 0. Therefore, D has a zero eigenvalue and owing 
to its symmetry a pair of real eigenvalues with opposite sign, namely ±|A;|. In view of ( ITT]) , this 
implies that Cl^ will also have a zero eigenvalue. Hence, its Jordan canonical form Jq2 must have 
one of the following two forms. 

Case 1 : J^^a = I A_ I , A± G C, (16) 




Case 2 : J^^ = A 1 , A G C. (17) 
\ A / 

These correspond to diagonalizable and non-diagonalizable cases, respectively. 



■^Given a linear operator acting in its representation in the Fourier basis, L, is related to its representation in 
the position basis, L, according to: LF{k) = (27r)~'^/^ e~^^'^ LF{x) dx^ . Here F is an arbitrary complex vector 

field having a Fourier transform F defined by F{k) := (27r)^'^/^ jj^g e^^'^'^ F{x) dx^. Throughout this article we use 
the same symbol for the representations of an operator in both Fourier and position bases to simplify the notation. 
This convention does not lead to any confusion, for we exclusively use the Fourier representation of the relevant 
operators. 



3 



Let S be an invertible matrix fulfilling Cl"^ = S ^Jj^zS*. Then, in view of (IMI) and (fTSjl . 



cos Wo 



tn) = s~'lj2 



.l)n 



vri=0 



{2n)\ 



[UJ, 



-V. 



n=0 



(19) 



In order to simplify these equations we consider Cases 1 and 2 separately. 
Case 1 is diagonalizable) : In this case, 



/ 

(Jn^r = 
V XI 

Inserting this relation in ( |T8l) and ( |T9|) yields 

cos(u;otf2) = S" ^ j cos{^,/Xl ujot) 
\ 

/ 

(l-^ sin(a;ot^) = S'^ 



for all n G Z"* 



sm{^,/\I Uot) 



\ 







s, 

COs(a/A^iW0^) / 








(20) 



S (21) 



Note that the time-harmonic solutions of the wave equation (fT3|) are the eigenvectors of Cl"^ 
[21 [6] . Hence, they are given by 



^±(A^)e*('=-''-V^±"°*) +^±(A?)e'('=-"^+V^±""*) S^^Ca, a = 1,2, (22) 



where and are complex-valued functions, and 



Clearly, for a = 1 and a = 2, S'^Ca is the second and the third column of the matrix 5*"^, 
respectively. Furthermore, the matrix S and the eigenvalues X± that enter (1221) are, in general, 
functions of k. 

Case 2 {Cl^ is non-diagonalizable): In this case. 















A" 


riA""^ 








A" 



for all n G Z"* 



(23) 
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and we find 



/ 1 



cos{uotCl) = S 



-1 











cosi^AuJot) -^sm{^AuJot) 
cos(a/Awo^) J 

/ 



(24) 



^sin(v^u;o^) ^ ujQtcos{y/XuJot) — ^sm{y/Xujot) ^ S. (25) 



sin(-\/A Wot) 



Moreover, the time-harmonic solutions of the wave equation f[T51) take the form 



E±(f,t) 



?f(A;)e 



i(A;-r'+\/A <^ot) 



(26) 



where ^ and $f are complex-valued functions. 



We close this section by revealing an implicit symmetry of the operator Q'^. Here we will 
momentarily consider the possibility that permittivity and permeability tensors are not constant. 

As we have seen, Q"^ admits a single zero eigenvalue. The presence of this zero eigenvalue implies 
the existence of parts of £ and fl that have no influence on the propagating electric field. In 
order to characterize these, we expand £ and /l as follows: 



1 «■/ -H-/ «■/ •(-!>/ 

£ = £Hq + £Hi + £Ao + 

-f^, •(-!>/ •<->/ 

= /^Ho+ /^Ao+ /^Ai- 



(27) 
(28) 



Here the subscripts H and ^4 denote the "Hermitian" and "anti-Hermitian" parts of the tensors 
respectively. 



and stands for the 3x3 zero matrix. 

Because the null space of is spanned by k, according to fl29|) we have 



(29) 



s'hi = cie(k) kk\ = il^e(k) kk\ = ci^(k) kk\ jlj^^ = i(3^(k) kk\ 

where a^, /?£, a^, are real-valued functions of k, and 

kl kik2 kik^ 



(30) 



kk^ 



For the cases that 



kik2 kj k2ks 
kik-i k2k3 kl 



■Ao - f^Ao 







(31) 



is pseudo-Hermitian. If in addition £ho is positive definite, then fi^ is quasi-Hermitian. Note 
that these hold irrespectively of the value of E//^, e^^, /i^^, and /i^^. 

The above observation seems to be relevant only for the material with spatial dispersion [TO] . 
For the non-dispersive material that we consider in this article, £ and /i are independent of k and 
the above gauge freedom is effectively frozen (the scalar functions a^, /?£, a^, are identically zero). 
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3 Applications: Plane- Wave and Time- Harmonic Solutions 



In this section we will examine the application of our general results in the study of three concrete 
models. This also allows for extracting the physical meaning of pseudo- Hermit icity and quasi- 
Hermiticity of the Q^. 

Example 1: Consider the propagation of a plane wave given by the initial conditions 

cosip 

Eoix,y,z) = S I sinv9 I e^'=^^ Eo{x,y,z) = 0, (32) 
in a uniaxial medium with permittivity and permeability tensors [9] : 




jji + i/3 

£ = eo \ ~ia ei + i^ | , = fJ^o \ -iP /ii + «7m | , (33) 

/i3 

where S & C and ip G [0, 27r) determine the amplitude and polarization of the initial wave, 
Eq and fiQ are the permittivity and permeability of the vacuum, so that eo/^o = ^/c^, and 
El, 62, /ii, /i25 7e5 7^*5 a, /3 G M describe the electric and magnetic properties of the medium^ 

A simple consequence of fl32|) is that in the calculation of in the Fourier basis we can set 

ki = k2 = 0. (34) 

This leads to an enormous simplification of the Jordan decomposition of Cl'^. In particular, 
Cl"^ is diagonalizable, uq = c|/c3|, and inserting (fTOl) and (l33l) in ( fTTl) and calculating the 
eigenvalues of the resulting expression for Cl^ we find 

A± = [{si ±a + Z7e)(/ii ±(3 + z7m)]^'- (35) 

Furthermore, the similarity transformation that diagonalizes is independent of the param- 
eters of the system, for we have 



S = \ i 10, = - 1 1 . (36) 





Substituting f p5|) and f l36|) in fl20|) and using the resulting relation and (|32!) in f|T3|) . we obtain 
the following expression for the propagating electric field 

E{f, t) = E{z, t) = EfiEit) e'^^\ (37) 

where 



^ / e-'^P cos( VA_ ujQt) + e*^ cos( VA+ wq^) 

V 



*In light of ([2]), Eo{x, y,z) = follows from Bo{x, y, z) — 0. 
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The expression for the magnetic field can be obtained using (jl]). Setting B{f,0) 
in view of ([3]) is consistent with the second equation in f p2|) . and using fl37|) and I 



0, which 
), we find 

(39) 



where 



nsit) := - 



e-"f 8111(^X7 wot) e''P sm(^A+ ujot) 







(40) 



It is remarkable that for the cases where is a pseudo-Hermitian operator but not quasi- 
Hermitian, so that A+ = Al, both the vectors nE(t) and ^^(t) stay real. Therefore, for 
the model we consider, the propagating wave does not undergo a phase shifljfl; the pseudo- 
Hermiticity of fl"^ implies phase conservation. In light of (15^ . A+ = \*_ and Q"^ is pseudo- 
Hermitian (but not quasi-Hermitian) if and only if /xia 7^ and 



ei/3 + i2ia = 0, £17^ + 



(41) 



In view of fl35|l . these relations imply 

X± = -[{el + ^!-a'T^a^e^'. (42) 
/il 

Next, we derive the necessary and sufficient conditions for the quasi-Hermiticity of fi^ (equiv- 
alently Cl"^). Because Cl^ is diagonalizable, the latter is equivalent to the reality of the eigen- 
values A-t of According to (135|) . A-t are real and is quasi-Hermitian if and only if 
either 

7. = 7m = 0, (43) 

or the following two conditions hold: 

eiP - i2ia = 0, £17^ + /ii7e = 0. (44) 

Condition P3l) corresponds to lossless media where £ and fl are Hermitian. In this case fi^ 
is known to be quasi-Hermitian [2] . Conditions (1441) do not require Hermiticity of either of £ 
and /i, yet they imply quasi-Hermiticity of Q"^. This is quite remarkable, because whenever 
dSD holds, 

A± = -[(ei±«)^ + 7a"'- (45) 
/il 

Therefore, X± are real and positive. According to ( 1371) - (HOl) . this means that the amplitude 
of both the electric and magnetic fields are bounded biperiodic functions of time; quasi- 
Hermiticity of confines the amplitude of the propagating electric and magnetic fields to a 
bounded region in the x-y plane that does not shrink in time. 

^According to p7|) . E(r,t) is given by a real veetor-valued funetion of x and t, namely fiE{x,t) times a eomplex 
sealar function of x, i.e., Ee^^^^ . Hence we can associate a total phase factor to E{r,t), namely £e^''^^/\£\ which 
does not change in time. In view of (|39|) . the same is true for B{f,t). 
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In view of (HT]) and fl44l) . a necessary condition for the pseudo-Hermiticity of Q"^ is the second 
equation appearing in both ( HTl) and (jHj). For the usual material where /xi and ei are positive, 
this equation implies that 7^ and 7^ have opposite sign. Therefore, to maintain the pseudo- 
Hermiticity of the material must simultaneously display loss with respect to the electric 
properties and gain with respect to the magnetic properties or vice versa. 

As seen from (138!) and (140!) . the propagating plane wave determined by the initial conditions 
( l32l) is not a time-harmonic solution of the Maxwell equations. It is rather a superposition of 
two time-harmonic plane-wave solutions. In general, the time-harmonic right-going, plane- 
wave solutions propagating along the 2;-axis has the following form. 

where A^^^ ^J^'^^ E C. These solutions are clearly circularly polarized. For the cases that Q,^ 
is quasi- Hermitian, i.e., either f H3|) or f H^ holds, fH6|) are periodic solutions with a constant 



amplitude. For the cases that is pseudo-Hermitian but not quasi-Hermitian, either E^_l^^ 
or E^^ is an exponentially decaying solution while the other is an exponentially growing 
solution. 

Example 2: Consider the propagation of electromagnetic waves in a medium with complex 
symmetric permittivity and permeability tensors |llj : 

£ = EqA, jl = /ioA, A := s b () , (47) 

V U [) c / 

where a, b, c,0, [),u G C. In this case, VL^ is a diagonalizable operator with a single nonzero 
eigenvalue, namely 

^± — ^0 •— • (4o) 

|A;|2[abc + 2q\)u - (af)2 + bu^ + cq^)] 

A direct consequence of this is the fact that pseudo-Hermiticity of VL^ coincides with its 
quasi-Hermiticity. The latter is achieved by requiring that Aq be real. 

We can easily diagonalize Vt^ by setting 



^3 afci+gfc2+ufc3 afei+gfc2+ufc3 



v 



I 1 I- (49) 

1 1 



In the following we derive the expression for the time-harmonic plane-wave solutions of 
Maxwell's equations for this case. In view of fl22|) and fH9l) . these are linear combinations of 

-(uki + i)k2 + cks) 

E^^\f,t) = 4^1 W^(l+)gi(fc-r-v/A^;^0i) +^{l-)gi(fc-r+v^^0i))^ (59) 




aki + + uk 



3 



. / -idh + bk2 + {)A;3) 



\k\ 
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where A^^"^^ and A^"^"^^ are possibly /c-dependent complex coefficients. 

Again whenever is quasi-Hermitian one obtains periodic time-harmonic solution. Other- 
wise, depending on the imaginary part of V^o one may have an exponentially decaying or 
growing solution. The latter is clearly dependent on the magnitude and direction of k. For 
example, for the special case of a plane wave propagating along the positive 2;-axis, we have 

= 1^ ^ e'('=^""^"«*\ 3^\z,t) = A^^+^ ( u ^ e'(^'^"^^""*\ (52) 

where uq = ck^ and 




ab + - — + 



c \ c 



(53) 



For the special case, a = b = (1 + u^)/c, q = u^/c, f) = u, that is considered in [TT], we 
have Ao = c^. Therefore, Q"^ is quasi-Hermitian and the plane wave solutions fl52l) do not 
decay in time provided that c is real. They are exponentially decaying (growing) solutions, 
for Im(c) > (Im(c) < 0). 

Example 3: Consider a time-harmonic plane wave propagating along the positive z-axis in 
a medium with complex symmetric permittivity and permeability tensors: 



£ =eo\ g f + iQ \ , }1 = l^ol, (54) 



where f and g are nonzero (possibly) complex parameters. In this case, ki = k2 = 0, uq = cks, 
and we can easily show that is a non-diagonalizable operator, i.e., it is an example of Case 2 
of Section 2. Moreover, we have 









M 






1 







\ 










A = f, S = 1 , S"^ = 1 I . (55) 
V g ^0 / \ 

According to fl26l) the propagating electric field is given by 

E{z,t) = A+\ % I e^(^'3--^'^ot)^ 




where G C. Again fi^ is pseudo-Hermitian provided that f is real. In this case, fl56p 
is periodic in time. Otherwise, its amplitude is an exponentially decreasing or increasing 
function of time. These correspond to Im(v^) > and Im(-\/f) < 0, respectively. 

Another peculiarity of the model considered here is that the expression fl56l) for the time- 
harmonic plane-wave solution does not involve g. We only require that g takes a nonzero 
value. Furthermore, ( 156|) describes a right-going clockwise circularly polarized field. The 
non-diagonalizability of ^ implies that there is no other time-harmonic plane wave solutions 
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propagating along the positive 2-axis that is hnearly independent of (156|) . We may take this 
as an indication of the "left-handedness" of the medium. This observation motivates the 
solution of the Maxwell's equation associated with the following initial conditions. 



E,{r) = E 



Bnir] 



(57) 



The result reads 



E{f,t) = E{z,t)=£ 



( cos(v^a;ot) + ^ sin(v^a;ot) \ 



% sin(v/f a;ot) 



ik^z 



B{f,t) = B{z,t)=c~^£ 



( 



\ 



sm 



(5^ 



(59) 



where we have made use of (ITSl) . ffT7|) . fl2^ . fl2S]) . and (jl]). The appearance of the factors 
Wot and bj\t^ as coefficients of the period functions of t in (158|) and (!59|) is a clear indication 
that these solutions are not time-harmonicj^ Indeed, even for the cases that f is real, these 
solutions are unbounded functions of time. This in turn means that one cannot maintain 
right-handed circularly polarized initial waves of the form f l57p in such a material|l| 



4 Concluding Remarks 

The information on the propagation of electromagnetic waves in a stationary, possibly inhomoge- 
neous, and anisotropic media is encoded in a generally matrix- valued differential operator In 
the absence of loss and gain, this operator is quasi-Hermitian and one can use the properties of 
quasi-Hermitian operator to obtain the solution of the wave equation. In the presence of gain or 
loss, there is no guarantee that this operator is diagonalizable or has real eigenvalues. 

In this article we examined the Jordan decomposition of this operator for a class of anisotropic 
active media and showed how the spectral properties of ^ might be used in the description of the 
propagating waves. In particular, we examined the time-harmonic and plane-wave solutions for 
three class of toy models, derived explicit conditions on the permittivity and permeability tensors 
that would render pseudo-Hermitian or quasi-Hermitian, and offered a physical interpretation 
of the pseudo-Hermiticity and quasi-Hermiticity of ^ in terms of the behavior of the propagating 
plane-waves solutions. 

An interesting observation we made was that ^ might actually be non-diagonalizable. We 
constructed an explicit model with this property and demonstrated a surprising feature of this 
model that could be interpreted as its left-handedness. Whether such exotic material can exist in 
nature or be manufactured is a subject of a separate investigation. 

Our results may be extended in at least the following two main directions. Firstly, one may 
try to generalize the method to non-homogeneous isotropic media. This can be achieved using 

^Note that g ^ 0. 

^There is not reason to believe that such exotic material can exist in real life. 
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the methods of Fourier analysis. Secondly, one might attempt to include the effects of dispersion 
following the ideas presented in [6]. 
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